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18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting # 7: The Fundamental Solution and Green Functions 



1. The Fundamental Solution for A in 

Here is a situation that often arises in physics. We are given a function f(x) on M" representing 
the spatial density of some kind of quantity, and we want to solve the following equation: 

(1.0.1) Au{x) = f{x), x = eM". 

Furthermore, we often want to impose the following decay condition as |x| ^ oo : 

(1.0.2) \u{x)\^0. 

For technical reasons, we will need a different condition in the case n = 2. A good physical example 
is the theory of electrostatics, in which u{x) is the electric potentiaQ and f{x) is the charge density. 
f{x) could be e.g a compactly supported function modeling the charge density of a charged star, 
and we might want to know how the potential behaves far away from the star (i.e. as \x\ oo). 



Roughly speaking, the decay conditions (1.0.2) are physically motivated by the fact that the star 



should not have a large effect on far-away locations. 



As we will soon see, the PDF (1.0.1) has a unique solution verifying (1.0.2) as long as f{x) is 
sufficiently differentiable and decays sufficiently rapidly as ^ oo. Much like in the case of the heat 
equation, we will be able to construct the solution using an object called the fundamental solution. 

Definition 1.0.1. The fundamental solution $ corresponding to the operator A is 

X X def f 77- In \x\ n = 2, 
^^•0-3) *(^)= --^ n>3 

where as usual |x| \/Er=i(^*)^ is the surface area of a unit ball in (e.g. uj^ = 47r). 

Remark 1.0.1. Some people prefer to define their $ to be the negative of our $. 

Essentially, our goal in this section is to show that A$(a;) = S{x), where 6 is the delta distribution. 



Let's assume that this holds for now. We then claim that the solution to ( 1.0.1 ) is u{x) = f * $(a;) = 
/Rn/(y)$(2; - y)dP-y. This can be heuristically justified by the following heuristic computations: 
Ax(/=^$) = />*-A,$ = />.(5 = /(x). 

Let's now make rigorous sense of this. We first show that away from the origin, the fundamental 
solution verifies Laplace's equation. 

Lemma 1.0.1. If x i^O, then A<l>(x) = 0. 



^Recall that the the force F associated to u is F = -Vu. 
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Proof. Let's do the proof in the case n = 3. Note that ^(x) = $(r) (r \x\) is spherically symmetric. 
Thus, using the fact that A = + ^dr when r > for spherically symmetric functions, we have that 
A$ = (92$ + ^5,$ = ^ + ^ = 0. ' □ 

We are now ready to state and prove a rigorous version of the aforementioned heuristic results. 

Theorem 1.1 (Solution to Poisson's equation in M"). Let f{x) e C|^(]R") (i.e., f{x) is a 
smooth, compactly supported function on W^). Then for n >3, the Laplace equation Au{x) = f{x) 
has a unique smooth solution u{x) that tends to as \x\ oo. For n = 2, the solution is unique 
under the assumptions ^^j^ as \x\ ^ oo and |V'u(x)| as \x\ oo. Furthermore, these unique 
solutions can be represented as 

(1.0.4) «(x)-($./)(x)-| __L/^J^|n-2;(^_^)rfn^, ^ > 3. 

Furthermore, there exist constants Cn > such that the following decay estimate holds for the 
solution as \x\ oo : 

, / N , f C2 In Ixl n = 2, 
(1-0.5) Hx)\<\ J^' ' ^>3' 

Remark 1.0.2. As we alluded to above. Theorem |1.1| shows that A$(x) = S{x), where 5 is the 
"delta distribution." For on the one hand, as we have previously discussed, we have that f = 6 * f. 
On the other hand, our proof of Theorem 1 1 . 1 1 below will show that / = Am = A($ * f) = (A$) * f. 
Thus, for any /, we have 5 * f = (A$) * /, and so A$ = 5. 

Proof. We consider only the case n = 3. Let's first show existence by checking that the function u 
defined in ( 1.0. 4[ ) solves the equation and has the desired properties. We first differentiate under 



the integral (we use one of our prior propositions to justify this) and use the fact that A^f^x -y) = 
Ayf{x - y) (you can easily verify this with the chain rule) to derive 

(1.0.6) AMx) = -— [ ^AJ(x-y)d^y = -— [ Ay f(x - y) d^y 



To show that the right-hand side of (1.0.6) is equal to /(x), we will split the integral into two 
pieces: a small ball centered at the origin, and it's complement. Thus, let -Be(O) denote the ball of 
radius e centered at 0. We then split 

(1.0.7) AMx) = -^ f ^,Ayf(x-y)d^y-— f ^,Ayf(x-y)d^y'^=I + IL 

^ ^ " ^ ^ 47r Jb.(o) \y\ ^ ^ An JBm \y\ 

We first show that / goes to as e ^ 0+. To this end, let 

(1.0.8) M sup \f{y)\ + |V/(z/)| + \AJ{y)\. 



Then using spherical coordinates (r, u) for the y variable, and recalling that d^y = du (where 
u e dBi{0) c ]R3 jg point on the unit sphere and duj = sinO d6 d(j)) in spherical coordinates, we have 
that 
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(1.0.9) |J|< f ^Ayf(x-y) d^y<M [ r du dr = 2e^7iM. 

Clearly, the right-hand side of (1.0.9) goes to as e 0+. 

We would now like to understand the second term on the right-hand side of (1.0.7). We claim 
that 



1.0.10) 



|/(x)-JJ|-0 



as e 0+. After we show this, we can combine (1.0.7), (1.0.9), and (1.0.10) and let e ^ 0+ to deduce 
that A^.n(a;) = f{x) as desired. 



To show (1.0.10), we will use integration by parts via Green's identity and simple estimates to 



control the boundary terms. Recall that Green's identity for two functions u, v is 



'1.0.11) / v{x)Au{x)-u{x)Av{x)d^x= / v\I f^u{a) - uV f^v{a) da. 

Jfl Jdfl 



Using (1.0.11) and Lemma 1.0.1, we compute that 



(1.0.12) 



r 1 1 . r 1 1 



Above, V^(o.-) is the outward unit radial derivative on the sphere dB^{Q). This corresponds to the 
"opposite" choice of normal that appears in the standard formulation of Green's identity for 5^(0), 
but we have compensated for this by carefully inserting minus signs on the right-hand side of 

jv(o-)[i| = ~[^) that \a\ = e on 95^(0), and that da = e^du on 95^(0), 



(1.0.12). Recalling also that V, 



it follows that 



i du. 



(1.0.13) - f ^Ayf(x-y)d^y = - f eu ■ (V f)(x - eu) du + f fix-eu). 

^ ' jBm\y\ JdBm JaBm ' 

Using (1.0.8), it follows that the first integral on the right-hand side of ( 1.0. 13[ ) is bounded by 
47rMe, and thus goes to as e -> 0+. Furthermore, since / is continuous and since /^^^(o-) Idu = An, 
it follows that the second integral converges to ATTf{x) as e ^ 0+. We have thus proved (1.0.10) for 
3. 



n 



To estimate \u(x)\ as |a;| oo, we assume that f{x) vanishes outside of the ball Bji{0). It suffices 
to estimate right-hand side of (1.0.4) when > 2R. We first note the inequality j^^^ < which 



holds for \y\ < R and \x\ > 2R. Using this inequality and (1.0.8), we can estimate right-hand side of 
(plbby 



;i.o.i4) 



1 

An 



1 



IBr{o) \x-y 



f{y)d'y 



M r 

2n\x\ Jbr{o) 



Id^y 



2R^M 
3\x\ ' 



and we have shown (1.0.5) in the case n = 3. 
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To prove uniqueness, we will make use of Corollary 4.0.4 which we will prove later. Now if u,v 
are two solutions with the assumed decay conditions at oo, then using the usual strategy, we note 
that w = u - V is a solution to the Laplace equation 



(1.0.15) Aw = 

that verifies -> as \x\ oo. In particular, w is a bounded harmonic function on M^. We will 



show in Corollary 4.0.4 below that w{x) must be a constant function. Furthermore, the constant 



must be since -> as \x\ oo. 



□ 



2. Green functions for domains Q 



Our goal in this section is to derive an analog of Theorem 1.1 on the interior of domains n c 
Specifically, we will study the boundary value Poisson problem 



(2.0.16) Au{x) = f(x), xeOcM'^, 

u(x) = g(x), X € dQ. 

Theorem 2.1 (Basic existence theorem). Let g be a bounded Lipschitz domain, and let g € 



C{dn). Then the PDF (2.0.16) has a umque solution u e C^{n) n C(fi). 



Proof. This proof is a bit beyond this course. □ 

Definition 2.0.2. Let Q c M" be a domain. A Green function in Q is defined to be a function of 
{x,y) eQ xQ verifying the following conditions for each fixed x eQ : 



(2.0.17) 
(2.0.18) 



AyG{x,y) = 6{x), 
G(x, cr) = 0, 



yen 
cr e d^l. 



Proposition 2.0.2. Let $ be the fundamental solution (1.0.3) for A in M", and let e M"- be a 



domain. Then the Green function G{x,y) for can be decomposed as 



(2.0.19) G{x,y) = <^{x-y)-^x,y), 

where for each x eQ, (f)(x,y) solves the Dirichlet problem 



(2.0.20) 
(2.0.21) 



Ay(f){x,y) = 0, 

(f)(x, cr) = $(2; - cr). 



y e n, 



Proof. As we have previously discussed, A$ = S. Also using (2.0.20), we compute that 



(2.0.22) Ay{<^{x-y)- y)) = A,<l>(x -y)- A,0(x, y) = 6(x - y). 

Therefore, $(x - y) - (f){x,y) verifies equation ( |2.0.17[ ). 



Furthermore, using (2.0.21), we have that ^{x - a) - (j){x,cr) = whenever a e dQ.. Thus, $(x 
y) - (j){x,y) also verifies the boundary condition (2.0.18). 



□ 
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The following technical proposition will play later in this section when we derive representation 



formulas for solutions to (2.0.16) in terms of Green functions. 



Proposition 2.0.3 (Representation formula for u). Let $ he the fundamental solution (1.0.3) 
for A in M", and let c he a domain. Assume that u e C^(f2). Then for every x £ fi, we have 
the following representation formula for u(x) : 

(2.0.23) u(x) = -y)Ayu(y)d'^y - / $(x - cr)Vjv(^)'u(cr) (ia + / ■u((t)V^(^)$(x - cr) rfd . 



single layer potential 
Proof. We'll do the proof for n = 3, in which case 



douhle layer potential 



We will also make use of Green's 



identity (1.0.11). Let B^{x) be a ball of radius e centered at x, and let fi^ '^=^ Q,\B^{x). Note that 



dile = dilu -dB^{x). Using (1.0.11), we compute that 



-Au{y)d:'y= / -V^u{a)-u{a)vJ. :)da 

fij \x - y\ Jan, \x - a\ ^\x - 

1 r" 1 

= / 1 TVi^u(a)da- / u(a)V fyi -, rlrfa 

Jdn\x-a\ Jan V|a;-(j|/ 

- / 7V<ru(a)da+ / u(a)S/^ri-, Ada. 

JdB,{x)\x-a\ JdB,{x) V|x-(t|>' 

In the last two integrals above, iV(cr) denotes the radially outward unit normal to the boundary of 
the ball B^{x). This corresponds to the "opposite" choice of normal that appears in the standard 
formulation of Green's identity, but we have compensated by adjusting the signs in front of the 
integrals. 



Let's symbolically write (2.0.24) as 
(2.0.25) L = R1 + R2 + R3 + i?4. 

Our goal is to show that as e j 0, the following limits are achieved: 

• L ^ -47r - y)Ayu{y) d^y 

• Rl An X single layer potential 

• R2 -An X double layer potential 

• i?3 ->0 

• RA^ -Anu{x). 

Once we have calculated the above limits, ( [2.0.23 ) then follows from simple algebraic rearranging. 
We first address L. Let M = maXy^^ Au{y) . We then estimate 



(2.0.26) 



[ ^^Au{y)d^y- [ -^Au{y)d^y < [ -^\Au{y)\d^y 
Jn\x-y\ Jn^\x-y\ JB^{x)\x-y\ 



<M 



d^y as e i 0. 



Ib4x) \x - y\ 

This shows that L converges to j^^Au(y) d^y as e i 0. 

The limits for Rl and R2 are obvious since these terms do not depend on e. 

We now address R3. To this end. Let M' = max^^j^ |Vm(?/)|. We then estimate _R3 by 
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(2.0.27) 

\R3\ , I 



\x - a\ 



V^M(a) 



da 



< [ -M'da= Ane^ xe^^M' -> as e i 0. 

JdB^(x) e ' — V — ' 

surface area of dB^ix) 



We now address i?4. Using spherical coordinates {r,9,(j)) e [0, oo) x [0,7r) x [0,27r) centered at 
X, we have that da - r"^ sva.6 dO dcf). Therefore, fg^^f^x) \x-a\^ ~ l4,e[0 2TT] leeio-K] ^dOdcp = in. We now 
estimate 



(2.0.28) 



— R4:-\-u(x)] = u(x) + — / u(a)S7Mf„J-, Ada 

r 1 

Jf 1 
(i 1^)^'^ 
dBAx) \ \x-aY' 



1 

47r 



1 

< — max \u 

47r aedB^x) 

< max \u{x) - u{a)\ ^ as e i 0. 

a£dB^{x) 



This shows that R4 -4ttu(x) as e i 0. 



□ 



Theorem 2.2 (Representation formula for solutions to the boundary value Poisson 
equation). The solution u to (2.0.16) can be represented as 



(2.0.29) 



u 



(x) = - f f(y)G{x,y)d^y- f gia) Vj^G(x,a) da. 
Jn Jan ■. . 



Poisson kernel 



Proof. Applying Proposition |2.0.3[ we have that 

(2.0.30) u{x) = - f ^{x-y)f{y)d''y+ f <^{x - a)S/ ff,„Ma) da - f g{a)S/ sj,„Mx - a) da. 

G(x,y) = ^{x-y)-(f){x,y) 



Recall also that 
(2.0.31) 



(2.0.32) 



G(x, a) = when a € dVl. 



Applying the Green identity (1.0.11) to the functions m(?/) and (f)(x , y) , and recalling that Ay(j)(x,y) 



for each fixed x £ fi, we have that 



Au{y) 



^(x-a) 



(2.0.33) 0= [ <j){x,y)f{y) d^y- [ 0(x, a) V^M(a) rfa + [ g{a) V j^(p{x,a) da. 

Adding ( |2.0.30D and ( |2.0.33D , and using ( |2.0.3ip , we deduce the formula ( |2.0.29[ ). 
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□ 



3. Poisson's Formula 



Let's compute the Green function G{x,y) and Poisson kernel P(x, a) -V;vG(x, a) from (2.0.29) 
in the case that n = Br(0)c M3 

is a ball of radius R centered at the origin. We'll use a technique 
called the method of images that works for special domains. 

Warning 3.0.1. Brace yourself for a bunch of tedious computations that at the end of the day will 
lead to a very nice expression. 



The basic idea is to hope that (j){x,y) from (2.0.19), viewed as a potential that depends on 
y, is equal to the potential generated by some "imaginary charge" q placed at a point x* e Q^. 



To ensure that property (2.0.18) holds, q and x* have to be chosen so that along the boundary 
{l/€M3 I \y\ = R},^{x,y) 

~ iiT\x~y\ • ■'■^ ^ nutshell, we guess that 



(3.0.34) 



G{x,y) 



At(\x - y\ 4:Ti\x* - y\ ' 



and we try to solve for q and x* so that G{x,y) vanishes when \y\ = R. Thus, when \y\ = R, we must 
have 



(3.0.35) 
(3.0.36) 



1 



4:7r\x - y\ 47r|x* - y\ 



X* - y\ = q \x - y\ . 



(3.0.37) |xp -2x-y + R^ = \x-y\'^ = q'^\x* - = q'^{\x*\^ -2x* ■y + R^). 

Then performing simple algebra, we have 



(3.0.38) 



xT + R - q {R + \x\ ) = % • {x* - q^x). 



Now since the left-hand side of (3.0.38) does not depend on y, it must be the case that the second 



term on the right-hand side vanishes. This implies that x* = g^x, and also leads to the equation 



(3.0.39) g>P - q\R^ + |xp) + R^ = 0. 

Solving ( [3.0.39 ) for q, we finally have that 



(3.0.40) 

(3.0.41) 
Therefore, 



R 

q = r-r 
\x\ 

. R^ 



X" 



-X. 



\x\ 
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(3.0.42) 



(3.0.43) 



my) 



1 R 



where we took a limit as x in (|3.0.42|) to derive (3.0.43). 

R 



Next, using (3.0.34), we have 
(3.0.44) 



(3.0.45) 



(3.0.46) 



G{x,y) 



1 



1 



X ^0, 



1 1 



47r|y| R- 
VyG{x,y) 



X - y 1 R X* - y 



47r|x-?/p 47r |x| |x* - 



Now when a e 03^(0), (3.0.36) and (3.0.40) imply that 



(3.0.47) 



\X - cr 



\X - cr . 



\X\ 



Therefore, using (3.0.46) and (3.0.47), we compute that 



(3.0.48) 



VaG(x,a) 



X - a 



1 |xp X* - cr 



X - a 



1 |xp ~ ^ 



47r|x - 0"^ At: R^ \x - 47r|x - At: R^ \x - 

12 , 



-cr 



At:\x - a\ 



-ii-m. 

\ R2 ) 



Using (3.0.48) and the fact that N{(j) = ^, we deduce 



(3.0.49) 



^N(a)G{x, cr) =^ VaG{x, a) ■ N{a) 



i?2 - bp 1 



Remark 3.0.3. If the ball were centered at the point p e M'^ instead of the origin, then the formula 



(3.0.49) would be replaced with 



(3.0.50) 



R^ -\x - pp 1 
AnR |x-crp 



Theorem 3.1 (Poisson's formula). Let Bji{p) W' he a hall of radius R centered at p = 
{p^,p^,p^), and let x = {x^,x'^,x'^) denote a point in M^. Then the unique solution u € C'^(Bji{p)) n 
C(Br{p)) of the PDF 



(3.0.51) 



j Am = 0, X eQ, 

1 u{x) = f{x), X 6 

can he represented using the Poisson formula: 
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(3.0.52) 



R"^ -\x - 
4^ 



L 



da. 



dBii(p) \x - a\ 



Remark 3.0.4. In n dimensions, the formula (|3.0.52|) gets replaced with 
(3.0.53) 



u{x) 



B? - \x -pp 



f 

JdBR( 



UnR JdBnip) \X - a 

where as usual, Un is the surface area of the unit ball in 



da. 



Proof. The identity (3.0.52) follows immediately from Theorem 2.2 and (3.0.50). 

4. Harnack's inequality 



□ 



Theorem 4.1 (Harnack's inequality). Letu be harmonic and non-negative in the ballBn^O) c 
M". Then for any x e -Bij(O), we have that 



(4.0.54) 



Proof. We'll do the proof for n = 3. The basic idea is to combine the Poisson representation formula 



with simple inequalities and the mean value property. By Theorem 3.1, we have that 

i?2 - bp 



(4.0.55) 



u{x) 



AttR JdBa{Q) \x - Cr|' 



da. 



By the triangle inequality, for a e 95^(0) (i.e. \a\ - R), we have that |x| - i? < |x - (t| < \x\ + R. 
Applying the first inequality to ( 4.0.55[ ), and using the non-negativity of /, we deduce that 



(4.0.56) 



u{x) < 



R + 


X 


1 


i?2- 


X 


^AttR 



l(0) 

Now recall that by the mean value property, we have that 



f f{<y)da. 



(4.0.57) 



u(0) = ^— [ fia)da. 



Thus, combining (4.0.56) and (4.0.57), we have that 



(4.0.58) 



u{x) < 



which implies one of the inequalities in (4.0.54). The other one can be proven similarly using the 
remaining triangle inequality. 

□ 

Corollary 4.0.4 (Liouville's theorem). Suppose that u e C2(]R"-) is harmonic on M". Suppose 
their exists a constant M such that u{x) > M for all x e M", or such that u{x) < M for all x e M". 
Then u is constant. 
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Proof. We first consider the case that u{x) > M. Let v'^= u + \M\. Observe that f > is harmonic 
and verifies the hypotheses of Theorem 4.1[ Thus, by (4.0.54), if a; e MJ^ and R is sufficiently large, 
we have that 



(4.0.59) 



i?'^-2(i?-|x|) 



i?'^-2(i?+|x|) 



- — \ , ' ' m(0) < u(x) < — — \ , ' : u(Q) 



Allowing cx) in ( |4.0.59 ), we conclude that v{x) = v{0). Thus, t> is a constant-valued function 
(and therefore u is too). 

To handle the case u(x) < M, we simply consider the function w(x) -u(x) + \M\ in place of 
v{x), and we argue as above. 

□ 
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